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Cosmic ray propagation is diffusive because of pitch angle scattering by waves. We demonstrate 
that if the high-amplitude magnetohydrodynamic turbulence with B / (B) ~ 1 is present on top of 
the mean field gradient, the diffusion becomes asymmetric. As an example, we consider the vertical 
transport of cosmic rays in our Galaxy propagating away from a point-like source. We solve this 
diffusion problem analytically using a one-dimensional Markov chain analysis. We obtained that the 
cosmic ray density markedly differs from the standard diffusion prediction and has a sizable effect 
on their distribution throughout the galaxy. The equation for the continuous limit is also derived, 
which shows limitations of the convection-diffusion equation. 


I. INTRODUCTION 

Propagation of charged particles, whose Larmor radii, 
rp, are much smaller than the field inhomogeneities is 
ballistic along the field lines, neglecting drifts. When 
small-scale turbulence is also present, it can induce pitch- 
angle diffusion, so the parallel propagation becomes dif¬ 
fusive, and the particles can also jump across fields lines 
through a distance rp. It is this type of propagation of 
low energy cosmic rays (CR) through the Galaxy that is 
believed to occur. Indeed, the characteristic field correla¬ 
tion scale. A, of the galactic magnetic field is of the order 
of a few parsecs mm- Hence, rp X for the CR below 
the knee, E < PeV, so they nearly follow field lines in 
a diffusive manner. The small-scale Alfvenic turbulence 
responsible for the pitch-angle diffusion is believed to be 
self-generated by the streaming cosmic rays [5]. 

The presence of strong magnetohydrodynamic (MHD) 
turbulence in the interstellar medium (ISM) [H |2] also 
results in chaotic distribution of the filed lines in the 
galaxy. Such a chaotic field line topology leads to an 
additional, very efficient three-dimensional (3D) diffu¬ 
sion of CR throughout the entire ISM m HI- More¬ 
over, the MHD (likely Alfvenic) turbulence in the ISM 
is known to be of high amplitude at the outer scale, that 
is B/{B) > I. Such high-amplitude magnetic field fluc¬ 
tuations affect diffusion via mirroring and transient trap¬ 
ping effects as well In addition, trapping in MHD 

turbulence is intermittent and transient because large- 
amplitude, quasi-coherent Alfvenic wave-forms (“mag¬ 
netic traps” or “magnetic bottles”) are not static but 
exist for a certain life-time — the Alfvenic time. 

Finally, the mean field in the Milky Way (and all other 
galaxies too) is non-uniform. There is a net gradient of 
B toward the galactic mid-plane and the center of the 
galaxy. The exponential scale-height, H, of the galac¬ 
tic disk region dominated by magnetic fields is about a 
kiloparsec [2]. Interestingly, particle diffusion with mag¬ 
netic trapping in the presence of the net field gradient 
has, to our knowledge, never been addressed before, ex¬ 


cept for our earlier paper, where we argued for such a 
possibility [8]. In this paper we demonstrate that the 
mean field gradient modifies diffusion drastically, so that 
it becomes asymmetric in which transition probabilities 
from a current particle’s position onto neighboring ones 
are unequal. We stress that asymmetric diffusion is not 
to be confused with the standard anisotropic diffusion, in 
which the probabilities are equal but the diffusion coeffi¬ 
cient can be inhomogeneous and anisotropic, in general. 


II. TRANSPORT IN A “MULTIPLE-MIRROR 
MACHINE” 

Here we first discuss the effect of magnetic trapping 
on the transport of particles using the example of a 
“multiple-mirror machine” illustrated in Fig. This 
model serves as a good toy model of statistically isotropic 
high-amplitude MHD turbulence in the ISM. Let us as¬ 
sume that the machine consists of a chain of V ^ I 
identical magnetic traps, each of length A, the low-field 
strength is Bq and the mirror field is > Bq-, hence 
the mirror ratio is i? = Bm/BQ. Efficient confinement 
requires that ri, <C A, so we assume it is satisfied. The 
angular size of the loss cones, 6, is sin^ 9 = 1/R. 

Next, we assume that there is pitch-angle diffusion of 
particles, so they can be trapped and de-trapped in the 
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FIG. 1. A set of N equivalent symmetric magnetic traps or 
“magnetic bottles”. The minimum of the magentic field. So, 
is in the center, the maximum field, Bm, is at the magnetic 
mirrors on both ends of the bottle. 
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magnetic bottles while they propagate from one trap to 
another. The most efficient regime is, of course, when 
the mean free path, Ipa^ due to the pitch angle diffusion 
(i.e., the length over which the pitch angle changes ap¬ 
preciably, by about 7r/2) is comparable to A. If otherwise 
Ipa ^ A, then the system is highly collisional, so the effect 
of trapping is suppressed. Alternatively, if Ipa S> A, the 
particles can move through several traps in one ballistic 
run, before they become trapped again. This effectively 
increases the trap size and diminished the effect of trap¬ 
ping too. For the sake of argument, we assume that the 
rate of pitch-angle diffusion is optimal: Ipa — A. For 
cosmic rays in the galaxy, the pitch-angle diffusion coef¬ 
ficient is energy dependent, so does the above relation. 
Here we also introduce the effective collision time-scale 
T Ipa/vth, where Vth ~ c is the thermal speed being of 
the order of the speed of light for CR, of course. 

Transport in such multi-mirror machines has first been 
studied elsewhere [3]. Obviously, such transport is dif¬ 
fusive: upon de-trapping from i-th trap, a particle will 
randomly go toward (z -I- l)-th or (i — l)-th one with 
50% chance. Next, we note that transport is the flux of 
streaming particles, i.e., the particles in the loss cones. 
Their fraction is about ~ 1/R, therefore the number of 
streaming particles is Ustream ~ no/R, where Uq is the to¬ 
tal particle number. The particle flux is ~ ristreamVth ~ 
noVth/R, that is, the effective plasma drift speed through 
the traps is ~ Vth/R- Finally, the effective diffusion co¬ 
efficient through the traps is suppressed by a factor of 
1/R^, namely 
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FIG. 2. The cartoon representing magnetic bottles in the 
presence of the mean field gradient. This is also a simplified 
one-dimensional model of the Galaxy. Here, the small waves 
in the trap centers stand for the low-amplitude Alfenic tur¬ 
bulence responsible for pitch-angle scattering of cosmic rays. 
The trap size, A, is the effective correlation length on the outer 
scale of Alfvenic turbulence in the Galaxy. 



where D ~ is the collisional diffusion coefficient in 
the absence of magnetic traps. We again mention that 
here we considered the most optimal case of Ipa — A, 
which can be written as A VthT/R <C A, where L is the 
size of the whole system. 

Using this result for CR propagation, we conclude that 
high-amplitude MHD turbulence present in the ISM sup¬ 
presses their diffusion coefficient by a factor of the order 
of 
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that is by about an order of magnitude if R/ (B) ~ 2 for 
example. 


III. ASYMMETRIC DIFFUSION 

Now, we turn to the case with non-vanishing mean field 
gradient. For this we modify the previous multi-mirror 
machine toy model by imposing a long-scale variation of 
the global field. The magnetic traps become asymmetric, 
as shown in Fig. This is the simplest, yet much more 


FIG. 3. The particle distribution function in a center of an 
asymmetric magnetic trap with Bq B 2 < B\, where Bq is 
the field at its center. Green wedges along U|| are the loss- 
cones with the corresponding opening angles 0i < 02, the 
blue wedges are redions of trapperd particles, red wedges at 
U|| < 0 are the “semi-trapped’ particles representing the “leak¬ 
ing out” population- they are reflected from the stronger field 
region on the left and then leak our from the trap through the 
right end. White arros represent the effect of the pitch angle 
diffusion which traps and de-traps particles inside the trap. 


realistic one-dimensional model of the ISM, in which the 
global field strength has a gradient toward the galactic 
plane (i.e., to the left, in Fig. [^. The characteristic 
gradient scale \s H ^ 1.5 kpc, which is the observed 
galactic magnetic scale-height. The trap sizes are of the 
order of the outer scale of MHD turbulence. They vary 
across the Galaxy but typically of the order of A ~ 1 — 
10 pc [T]. 

The most important difference with the previous model 
is the slight asymmetry of the mirrors of individual traps, 
which leads to the substantial difference in the particle 
distribution function, see Fig. An asymmetric mag¬ 
netic trap has mirror fields Bi > B 2 ^ Bq, where Bq 
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FIG. 4 . The Markov chain representing the toy model of 
one-dimensional asymmetric diffusion. The transition prob¬ 
abilities to the left, r, and to the right, g, are constants but 
generally r ^ g. The in-flux of particles, Jin, into the sys¬ 
tem occurs through the right boundary. Each cell represents 
a magnetic trap; there separation is of the order of the cell 
size. A, which plays a role of the effective mean free path. 


is the field at its center. Correspondingly, the loss cones 
differ too, 62 > Oi. The loss cones, shown in green, cor¬ 
respond to the streaming population of particles. The 
trapped population (shown in blue) is, in contrast, sym¬ 
metric. Finally, the red population are the particles mov¬ 
ing to the left. They will be reflected from the left mirror 
and, in a half of the bounce time, will leak out from the 
trap through the weaker right mirror. White arrows il¬ 
lustrate pitch angle diffusion. 

Obviously, the particle populations escaping through 
both ends are not equal and proportional to the volumes 
of the spherical cones of the with the opening angles, 6 * 1 ^ 2 , 
namely 
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where Ai?i = Bi — Bq and AB 2 = B 2 — Bq and we 
used that cos 0 = (1 — sin^ 0 )^/^ = (1 — BojB)^^'^ = 
(AB/B)^^^. 

To proceed further we introduce a mathematical model 
of a one-dimensional Markov chain [ini E]. It is shown 
in Fig. 1^ Here, every magnetic trap is represented by 
a site, where a particle stays for some collision time, t, 
before it de-traps and rapidly jumps left or right into a 


neighboring site. As for any diffusive process on a Markov 
chain, we assume for simplicity that such de-trapping/re- 
trapping events occur at regular intervals, t. The num¬ 
ber of sites is A = LjX with L being the system size. 
In order to account for the leakage from the system, we 
have to add two “absorbers” at both ends of the chain — 
the ^-states, often referred to as the “limbo”-states. Fi¬ 
nally, we have to have a source of particles. For example 
and without loss of generality, we allow for the in-flux of 
particles, Jin, through the right boundary. 

The transition probabilities at each time step are r 
and g, respectively. Standard diffusion is known to have 
equal probabilities to go left or right — this is the classi¬ 
cal mathematical problem, the “drunken sailor problem”. 
In asymmetric diffusion, in contrast, the probabilities of 
the forward and backward transitions are not equal. In 
general, these probabilities may be different for different 
sites. Here we assume that they are constant throughout 
the chain. The values of r and g are the fractions of the 
left- and right-streaming particles, respectively, that is 
they are r = Hi/(r2i -|- XI 2 ) and g — 02/(Hi -|- H 2 ). It is 
convenient to cast them into the form: 

r=l/2-e, g = l /2 + e, (5) 

where 
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and SB = Hi — H 2 . Since the mean field varies only 
slightly on the scale of a single trap, we can estimate 
that 
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SB - A|VH| - {X/L)Bo. (7) 

The factor X/L can be interpreted as the ratio of the 
turbulence outer scale to the characteristic mean field 
gradient scale. We can also estimate that AB 2 = {B 2 — 
Bo) ~ B and Bi/Bo — B 2 /B 0 ^ R, so that 
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Note the appearance of the \/R^ suppression factor akin 
the one obtained in the previous section. 

The dynamics of a Markov chain is determined by the 
“continuity equations” written for each site. They are 

dtni =-(ri+gi)ni+r 2 n 2 , (9) 

dtn2 =-{r2 + g2)n2 + gini + Tono, ( 10 ) 


diUj = -{vj + gj)nj + gj-iuj^i +rj+inj+i, ( 11 ) 

dtnpf = —{rpf + gN)nN + gN-inN-i + Jin, ( 12 ) 

where dt stands for the time derivative and rij is the 
probability of a particle to be at the j-th site. Note that 
the first and last equations describe the left and right 
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FIG. 5. The solution of the one-dimensional Markov chain 
model with constant coefficients g and r. Asymmetric diffu¬ 
sion corresponding to the asymmetry parameter a = gjr ^ 1 
is shown by red curves. For comparison, the linear density 
gradiant predicted by the standard diffusion, a = 1, is shown 
with the blue curve. Here we took N = 100 and normalized 
al the curves to the same value at fc = 50 to ease comparison. 


defines the yet unknown ni. 

The final result — the analytical solution of the set of 
Markov chain equations (||)-@ — is 

nk = Jin ( Q,Af+l _ 1 ) ’ 1 < fc < iV. (18) 

This result is presented in Fig. which shows the par¬ 
ticle densities as a function of coordinate (to be precise, 
of the trap number, k). The salient feature of asymmet¬ 
ric diffusion is the non-linear particle density gradient. 
Depending on the value of a, one sees exponential sup¬ 
pression/enhancement of the density. The linear gradient 
is recovered in the case of symmetric, i.e., standard, dif¬ 
fusion for which a = 1. 


IV. GALACTIC EXAMPLE 


boundaries of our system. For an ensemble of CR parti¬ 
cles, rij is also the relative number density of particles in 
the j-th trap. Hereafter we assume the steady state, for 
simplicity, so that all time derivatives vanish, dtru = 0. 

We solve the system of N coupled equations (|^ 
by calculating the generating function [101E] defined as 
follows 

N 

(13) 

j=o 


where ^ is a dummy variable having no physical meaning 
whatsoever. Then, Uk at some fc-th site is simply the 
fc-th derivative of F{£^) evaluated at ^ = 0, namely 




1 




(14) 




The generating function is obtained, in turn, from the 
Markov chain equations as follows. One multiplies, the 
first equation by the first power of the second equa¬ 
tion by and so on until the last equation being multi¬ 
plied by Then we sum up all the resulting equations. 
Upon straightforward manipulations to isolate the com¬ 


binations obtain 


- 3=0 


m = 


- ni) 

[(1-b a)^ - 1 - 


(15) 


where we introduced the asymmetry parameter a = g/r. 
Using Eq. (14), we obtain 


Uk — (1 -|- Q! -|- Q.^ I 


0-^ — 1 
^) ni = -- ni (16) 

a — I 


valid for 1 < k < N. Next, we apply the right boundary 
condition, namely that the {N -b l)-th site is a limbo 
state, so riAT+i is identically zero. Thus 


- 1 


(17) 


We now give an illustrative example pertinent to cos¬ 
mic rays: the asymmetric CR diffusion in the direction 
orthogonal to the galactic plane. Fig. illustrates the 
Markov chain for this system. The Galaxy is seen edge- 
on and the galactic plane is vertical and located between 
sites N/2 and N/2 + 1. The source of CR is a supernova 
located inside the galaxy at some site with the number 
kNS = lV/2 -b m, that is m-th site to the right from 
the plane. The CR particle flux, Jcr, is ejected in both 
directions. As we’ve mentioned earlier, the size of the 
system is set by the exponential scale height, L — 2H 
with F[ ^ 1.5 kpc and the size of magnetic traps be¬ 
ing comparable to the outer scale of ISM turbulence. 


system size 



FIG. 6. The Markov chain representing GR diffusion in the 
Galaxy. The Galaxy is seen edge-on with the galactic plane 
being vertical. The galactic north is to the right. The sys¬ 
tem size is set by the galactic magnetic field scale-height, so 
the strongest field is at the mid-plane. Thus, the g and r 
probabilities are unequal on either side and, by symmetry, 
g+ = r_ and g- = r+. The cosmic rays are injected by a 
supernova (SN) located at some m’th site to the right, i.e., at 
ksN = n/2 -b m site. 


Rn+i — ~Jin + 


a — 1 


ni = 0 
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Upon doing so, the transition term from the (j — l)-th 
site to the j-th results in the following substitution 


^ ri{x)g{x) 


d{ng) 1 (P{ng) 
dx 2! dx"^ 


1 d^jng) 

3! dx^ 

( 20 ) 


Similarly, we treat the (j + l)-th site. 

Eq. (191 readily becomes, assuming r and g being 
constants. 


FIG. 7. The solution of the one-dimensional Markov chain 
model of our Galaxy with N = 2000, and the cosmic ray 
source located 300 cells away from the Galactic mid-plane 
toward the galactic north (i.e., to the right), that is at fcsjv = 
1300. The standard diffusion prediction (a = 1) is shown 
with the dashed line. The orange, blue and magenta curves 
represent the cosmic ray densities for asymmetric diffusion 
with Q = 1 ± e with e = 10“®, 10“^, 10“^ respectively. 


A ~ 1 — 10 pc. Thus, we estimate the chain length to 
be ~ L/\ - 1000. 

The magnetic field is inhomogeneous: it decreases 
away from the galactic plane toward both ends. There¬ 
fore r and g probabilities differ on the left (—) and on 
the right (-I-) from the plane. By symmetry, g_|_ = r_ 
and g- = r+. 

This Markov chain can be solved analytically in the 
same way outlined in the previous section. Alternatively, 
one can use the generic solution, Eq. (18), to piece-wise 
construct the global solution by matching the solutions 
with a < 1 and a > 1 at the mid-plane and the SN lo¬ 
cation. The resulting equation is cumbersome but not 
instructive, so we just plot the solution in Eig. This 
solution has N = 1000, ksN = 300 and several a = 1 ± e 
with e = 0, 10“®, 10“^, 10“^. Which one does corre¬ 
spond to our Galaxy? Assuming that the asymmetric 
diffusion-induced anisotropy should not exceed the ob¬ 
served CR anisotropy, which is about 0.1%, we estimate 
e to be e < 10“^. This corresponds to the blue curve, 
which shows about 50% difference compared to the stan¬ 
dard diffusion picture. 


V. COMMENT ON A CONTINUUM LIMIT 

It is instructive to determine how a finite-step-size dif¬ 
fusion on a Markov chain can be extrapolated to a con¬ 
tinuum limit. To do this, we consider the generic Markov 
chain equation 

diij ■ 

= rij-igj-i + rij+iTj+i - Uj{gj + Tj) (19) 

and treat j as a continuous variable — a dimensionless 
distance, x, measured in the units of a diffusion mean- 
free-path, that is A. Time is also dimensionless here be¬ 
cause it is measured in the units of the collisional time- 
step, T. As long as 1, one can Taylor expand in small 
A_) = 1 ^ J and substitute j with a;, for convenience. 


dn 

1h 
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(r + g) d'^n 
2 dx'^ 


(r - g) d^n 
3! dx^ 


■ • ( 21 ) 


The first two terms here are the constituents of the fa¬ 
mous convection-diffusion equation 


dn 


dn „d'^n 


( 22 ) 


Thus, the convection-diffusion equation is the expansion 
of the Markov chain equation up to the second order. 
Hence, it holds only when the higher order derivative 
terms are small enough. The drift (flow) velocity and 
the diffusion coefficient are identified as follows 


= r-g, D = {r + g)/2. 


(23) 


We also note that the density profile, Eq. (18), becomes, 
in the continuum limit and large x: 


n{x) oc 

if a > 1, that is exponential, and 

n{x) cx 1-6^*'^“, 


(24) 


(25) 


if a < 1. 


VI. CONCLUSIONS 

We demonstrated that asymmetric diffusion of parti¬ 
cles can occur in high-amplitude MHD turbulence with 
non-zero mean magnetic field gradient. The predicted 
density profiles are markedly different from the standard 
diffusion paradigm. We argued that accounting for asym¬ 
metric diffusion of CR in our galaxy can change their pre¬ 
dicted density distribution by a factor of two. We also 
indicated that the familiar convection-diffusion equation 
of particle transport is just the first two terms of expan¬ 
sion of the more generic Markov chain equations. 

We stress that the asymmetric diffusion discussed here 
should not be mixed with anisotropic diffusion. The for¬ 
mer deals with different probabilities of a particle to jump 
left or right, whereas the standard (even anisotropic) dif¬ 
fusion has those probabilities equal, though they may be 
position- and orientation-depend. 

Second important point to stress is that asymmetric 
diffusion results in a non-zero particle flux through the 
system, even though the plasma and the magnetic fields 
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are completely static. This is in contrast with a conven- the recently observed “Fermi bubbles”. 
tional lore that the global (“bulk”) motion of CR in the 
galaxy may only be associated with global motions of the 

ISM, i.e., winds. This may particularly be important for ACKNOWLEDGMENTS 
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